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In water infrastructures design problems, small changes in 
their geometries lead to a major variation in the construction 
time and costs. Dams are such important water 
infrastructures, which have different types regarding their 
materials and their behavior to endure loads. In the current 
paper, invasive weed optimization (IWO) algorithm is 
employed to find the best shape of a concrete gravity dam 
(Tilari Dam, India). Stress and stability were considered as 
design constraints, based on the following models: Model I 
(M1): upstream dam face is inclined and Model II (M2): 
upstream dam face is vertical. Optimization using IWO for 
M1 showed 20% reduction in cross-sectional area as 
compared to prototype. Although results obtained using [TWO 
showed no changes in comparison with the algorithms in the 
literature (i.e., differential evolution, charged system search, 
colliding bodies optimization, and enhanced colliding bodies 
optimization), it converged faster. But results for M2 
revealed 26% reduction in cross-sectional area. 
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1. Introduction 


Dams are essential infrastructures, which are built all over the world for meeting various water 
demands, including flood control, water supply (urban, domestic, industrial etc.), electricity 
generation, recreational activities, navigation, groundwater recharge etc. 


It is estimated that the gravity dams are the first water barriers in the history of human lives. A 
gravity dam is a heavy structure, which is made of concrete or masonry materials across the river 
to increase the volume and height of water. In fact, gravity dams are among the most common 
types of concrete dams that have received special attention because of their simple design and 
their applications in different types of valleys. The stability of a concrete gravity dam entirely 
depends on its mass. Normally, the weight of a gravity dam suffices for stability against all 
design loads. Although gravity dams have been constructed in different shapes, they are 
generally made with roughly triangular cross-sections [1]. They had been built with masonry 
materials before the 1800s [2]. Nowadays, they are mostly constructed with concrete. 
Trapezoidal and rectangular profiles were used to build the first samples of gravity dams’ cross- 
sections. Although the recent dams’ shapes have emerged by the development of new materials 
and design techniques, which aim to find more optimal shapes by researches and civil engineers. 


Optimization is an interesting technique in hydraulic structures design, which aims to find the 
best solution by searching the design variables in the search space [3]. Many studies in water 
engineering have been performed using various intelligent techniques, including spillways [3-7], 
reservoirs [8—10], earth dams [11], evaporation [12]. 


Generally, the structural optimization problems can be divided into three categories:(I) size 
optimization, (II) shape optimization, and (II) topology optimization [13]. Optimal design of 
dams can be done using mathematical methods or intelligent techniques like nature-inspired 
algorithms. The nature-inspired algorithms or in general metaheuristics tend to be better than the 
traditional methods on challenging, real-world problems due to the following reasons: 


e As traditional algorithms are mostly local search and gradient-based, so there is no 
guarantee for finding the global optimum due to the existence of a large number of local 
solutions in real-world problems. Consequently, the final solution will often depend on the 
initial starting points. 

e As traditional algorithms normally employ some information like derivatives about the 
local objective, they tend to be problem-specific. 

e Traditional algorithms are not able to solve highly nonlinear, multimodal problems 
efficiently, and they struggle to cope with problems with discontinuity, especially when 
gradients are needed. 


Different techniques have been employed to optimize gravity dams. Salmasi [14] optimized a 
gravity dam section using the genetic algorithm. Khatibinia and Khosravi [15] solved shape 
optimization problem of a concrete gravity dam using an improved gravitational search 
algorithm. Deepika and Suribabu [16] used Differential Evolution (DE) algorithm in order to 
find the best optimal shape of a gravity dam. The best solution was compared with an analytical 
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model and the results showed about 20% reduction in concrete usage of dam. Kaveh and Zakian 
[13] optimized a concrete gravity dam section using Charged System Search (CSS), Colliding 
Bodies Optimization (CBO), and Enhanced Colliding Bodies Optimization (ECBO) algorithms. 
The results were compared to DEA results of Deepika and Suribabu [16]. All three used 
algorithms had superior results to DEA. Chiti et al. [17] optimized a gravity dam shape subjected 
to earthquake load based on reliability—based design optimization. Khatibinia et al. [18] used the 
hybrid of an improved gravitational search algorithm and the orthogonal crossover to optimum 
design of concrete gravity dams. Memarian and Shahbazi [19] used DE algorithm in 
optimization of some gravity dams’ prototypes under various constraints. Zhang et al. [20] 
studied shape optimization of high RCC gravity dams regarding hydraulic fracturing. Apart from 
the aforementioned works, numerous problems in civil engineering have been solved using 
intelligent methods or, in general, soft computing methods [21-24]. 


Considering the No Free Lunch theorem [25] in optimization, which logically proves that there is 
no optimization algorithm to solve all optimization problems, designers need to evaluate various 
algorithms on a specific problem to see if it is better than others or not. Therefore, in the present 
study, a nature-inspired algorithm (1.e., invasive weed optimization algorithm) is employed to 
solve shape optimization problem of a concrete gravity dam based on two models geometry, 
considering a real benchmark design problem (i.e., Tilari Dam in Maharashtra, India). The 
selected gravity dam was optimized using some evolutionary algorithms in the previous works, 
which their results are compared to the current findings. Besides, in the present research, design 
variables bounds are changed to find better solutions for shape optimization problem of the 
concrete gravity dam. 


2. Methodology 


Nature-inspired approach, as a branch of artificial intelligence (AI), was chosen in the present 
study because intelligent methods are far faster and more precise than traditional methods and 
have viable results in the previous works. Optimization algorithms have different parameters, 
which should be determined at first step of optimization. These parameters are calculated using 
sensitivity analysis. The mathematical model of the problem is built considering the major 
factors and all design parameters. This model includes objective function of the problem (i.e., the 
area of gravity dam cross-section), design variables, and constraints. Penalty function technique 
is employed to consider them into objective function. 


2.1. Invasive weed optimization 


Invasive Weed Optimization (IWO) algorithm is one of the nature-inspired algorithms, which 
inspired by colonizing weeds and was introduced by Mehrabian and Lucas [26]. Comparison of 
the results of the IWO with four types of Evolutionary Algorithms (EAs) such as Genetic 
Algorithms (GAs), Memetic Algorithms (MAs), Particle Swarm Optimization (PSO) and 
Shuffled Frog Leaping Algorithms (SFLA) showed superior performance and convergence rate 
etc. [26]. Efficiency of TWO in optimization has been proved in different studies in water 
engineering [27-30]. 
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The process of achieving the optimal solution in the [WO is as follows: 


I. Initializing a population 


The implementation of this algorithm begins with the distribution of a certain number of seeds 
(initial population) in the search space. 


Il. Reproduction 


Each seed grows according to its merits and produces new seeds. The number of seeds produced 
by each plant increases linearly from the lowest possible number of seeds to the highest possible 
number. 


I. Spatial dispersal 


In this section, the generated seeds are randomly dispersed in the multidimensional search space 
by the normal random distribution. Its average value is zero and its standard deviation varies at 
different stages. This step is similar to the random propagation of the seeds around the parent 
plant. At each step, the value of the standard deviation o corresponding to the random function is 
reduced from the initial value of Oiniziqi to the final value of ofinai. In the simulations, the nonlinear 
change expressed in Eq. (1) has shown a performance: 


« Ger. ter)” 


max 


(iter,,,.)" 


oO, 


iter 


(Oinitiat _ O Final ) + O yinal (1) 


In Eq. (1), ité7max represents the maximum number of iterations, oj, the standard deviation in the 
current time step, and the nonlinear modulation index. 


IV. Competitive exclusion 


If the plant produces no seed, it will become extinct and otherwise, it can spread throughout the 
world. Therefore, some competition is needed to limit the maximum number of plants. After 
several iterations, the number of plants will reach their maximum. At this stage, it is expected 
that the more competent plants will proliferate than the other plants. When the number of plants 
reached its maximum (Pyax), the process of removing plants begins with less fitness [26]. 


On the other hand, in the current study, [WO performance is compared to four algorithms in the 
literature i.e., Differential Evolution (DE), Charged System Search (CSS), Colliding Bodies 
Optimization (CBO), and Enhanced Colliding Bodies Optimization (ECBO). DE is one of the 
widely used metaheuristics, which can be categorized as a population based optimization 
algorithm. DE was developed by Storn and Price [31] and has been used in many real 
engineering complex problems by employing mutation, crossover, and selection operators' 
techniques. CSS is a population based algorithm, which was developed by Kaveh and Talatahari 
[32] and was inspired from Coulomb law from electrostatics and the Newtonian laws of classic 
mechanics. CBO algorithm is a metaheuristic, which was presented by Kaveh and Mahdavi [33]. 
Unlike the most of optimization algorithms, CBO works simply and does not depend on any 
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internal parameter. ECBO is an improved version of CBO, which was proposed by Kaveh and 
Ichi Ghazaan [34]. Colliding memory was employed in this algorithm in order to collect the best 
solutions. Techniques in harmony search algorithm were used to improve CBO. 


2.2. Gravity dam optimum design model 


Figure 1 shows the schematic of a gravity dam (plan and section). The purpose of shape 
optimizing of a structure is to find the most appropriate dimensions and shape so that it can 
withstand all loads and pressures. The loads in the gravity dam models are divided into two 
major categories including vertical and horizontal loads. The vertical loads include self-weight, 
uplift pressure force, silt pressure force, and seismic force. In addition, the vertical loads include 
water force, silt pressure force, wave pressure force, and seismic forces. Sliding and shear failure 
occur when the horizontal forces on each horizontal plane of a dam exceed its shear strength. 
Overturning of the dam and additional compressive stresses (and possibly tensile) can be 
prevented by selecting the appropriate cross-section. Normally, a gravity dam may be failed due 
to one or all of these reasons: 


1) sliding on a horizontal plane 
2) overturning on toe 
3) weakness in material (stress > allowable stress) 


Design constraints are normally related to safety consideration according to the nature of the 
engineering problems and design codes. The architectural and usability issues may be also 
considered as design constraints, but these issues are generally considered as the solution ranges 
of design variables to constrain the generation of possible optimum solutions. An optimization 
problem requires objective function(s) or cost function(s), which is widely related to the cost of 
the design, but safety, usability and architectural problems can be added into the formulation 
[35]. The design variables directly affect the objective function, which their values are unknown 
at the beginning of solving problems. Evolutionary algorithms based on the optimization process 
obtain their values. The calculated values of design variables must be within a desired range. The 
optimal solution is acceptable in case it contains all constraints and limitations in design 
problem. To apply constraints, penalty function is usually employed to consider them into an 
objective function [3,36,37]. 


The objective function in optimization of the gravity dam’s shape is formulated as follows: 


Minimize A,,,, =0.5x,x,+BH +0.5x 3x, (2) 
where Adam is the area of gravity dam cross-section (m’) and other parameters are shown in Fig. 
1. The fitness function is formulated as Eq. (3) to include penalty functions (constraints) into the 
objective function (Eq. (2)). 


Fitness Function = A,,,, + > f,¥ (3) 


i=1 
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where f, = penalty values; Y= a large integer number which makes the objective function 


unacceptable (in the case of penalty 4 0); and m = number of constraints. 
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Fig. 1. A gravity dam schematic including geometric parameters. 


Various loads affect gravity dams design. They can be categories into two major group: vertical 
and horizontal loads. The forces in the gravity dam design can be presented as follows: 


A. Vertical forces 


1) Self-weight: force x (liver arm about toe) 
1 1 
Was itis x pare (4) 
B 
W,=y7,BH «(.43) (5) 
1 2 
W, = ated aha x Ga (6) 
| 2 
P, ~ awn x [:48 +3s,] (7) 
3 
1 
Py =%,% (A -x,) x [x48 +32) (8) 


1 ee mh 
R= lem Me x | (9) 


A. Ferdowsi et al./ Journal of Soft Computing in Civil Engineering 4-3 (2020) 65-78 


2) Uplift pressure force: force = (liver arm about toe) 


1 ; 
U, =3% (x, +d,)(h-h) x [s+ +30%,-4,) 


1 1 
U, = 3% (x, +d, (h+2h) x [+8 Heed) 


1 h-h 2 
Us = 5M Xs +B -d, )( ; ) x Fo. B ‘| 


U,= 4, (05 +B -d,)h a fe +B -4,) 


3) Silt pressure force: force x (liver arm about toe) 
= 1 2 nh 
Fy, =~x0.925y, nh; x +, bB x= : 
2 3) 
4) Seismic force: force x (liver arm about toe) 
EV, =aW, x C +B+— 7, | 
EV,=aW, x [:,4 +) 
EV, =aW, « (2: | 
3 
EV,=a,P, x hs, B42s,] 
1 
EV;=a,F,, x [. +B +n) 
v,=aR, x(t] 
B. Horizontal forces 
1) Water force: force x (liver arm about toe) 


1 h 
P = hh x }— 
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(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


(22) 
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2) Silt pressure force: force x (liver arm about toe) 
I sec file 
Py = eal ho x 3 (23) 
3) Wave pressure force: force x (liver arm about toe) 
3 
By =2y,h,°x U tah (24) 
4) Seismic forces: force x (liver arm about toe) 
EH, =a,W, +) (25) 
EH, =a,W, x (=) (26) 
Xx 
EH, =a,W, x [ (27) 


P,, =0.726p,,h 
Pex =C,,4y7,h (28) 
M,, =0.299C a,7,,h? 


P= 0.126p,,,h 
Pay =C ny toh (29) 


M , =0.299C',a,7,h° 


In optimization and design of structures, values of some parameters are fixed. In fact, they have 
been chosen by the experience or using previous experiments, design, and experiences, codes 
and by designer judgment. These fixed parameters in the gravity dam problem of Tilari Dam are 
as follows: 


Dam height (# )= 38.55 m 

Maximum (upstream) water level (h )= 36.2 (m) 

Maximum (downstream) water level (h')= 3 (m) 

Silt deposit level (h, = 13 (m) 

Specific weight density of water (7, )= 9.81 (kKN/m3) 

Specific weight density of concrete (7, )= 2.47, 

Friction coefficient of (1 )= 0.75 

Permissible shear stress at foundation (4 )= 1200 (kPa) 
Permissible compressive strength of concrete (c,)= 3000 (kPa) 


S008 al. OY, Ne OES 
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10. Crest width (B )= 4.9 (m) 

11. Downstream face height (*,)= 33.35 (m) 

12. Fetch (f )= 10 (km) 

13. Wind velocity (V,, = 80 (km/h) 

14. Centre of drainage gallery from axis (d, )= 1 (m) 


Five variables are selected as design variables. These variables and their upper and lower bounds 
are represented in Eq. (30). Two models are considered in present paper. Normally, the upstream 
and downstream slopes (n and m) are considered between 0-0.2, and 0.6—0.8, respectively [2] . 
These parameters in model I (M1) were chosen as 0.1-0.2 and 0.6-0.9 according to studies of 
[13] and [16]. In model II (M2) the upper face of gravity dam is considered perpendicular (n=0). 
O.1<n <0.2 
0.6<m <0.9 
Design variables  0.8h <x, <0.95h 
0.05 <a, <0.2 
0.05 <a, <0.2 


(30) 


Stability (overturning and sliding), stress, and geometry constraints are applied in shape 
optimization in the current study. The stability and stress constraints are shown in Eqs. (31) and 
(32). The geometry constraints are applied to the problem using upper and lower bounds on 
design variables. 


M 
Overtuming — FOS = 2M q 21:5 


Oo 
ME, 


H 


Stability ; Sliding > FSS = 21 


F+gB 
Shear Friction Factor > SFF = fee >3 


Constraint s 
- (31) 
Pee 0S 6, 
Principal 
Heel 40, So, 
Stress 
Toe > Typ SO, 
Shear 
Heel > T,,y S09, 
in which 
On = yp sec” Pp - Oy — p_,,)tan” y 
Oy = Oy sec’ d, —(Py — Py ) tan’ g, 
(32) 


Typ = Kom _ (Px a oper )] tan ¢, 


Thu = lov _ (Pi + Pen )] tan d, 
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SR. be . SF . be 
= (ae 5 Py =%,h 4 ae - =— (+=), and Py =¥,h. 


1 1 1 1 


where o,, = 


3. Results and discussions 


Sensitivity analysis for choosing the [WO parameters was conducted and its results are shown in 
Table 1. In fact, sensitivity analysis is necessary to gain the best value of the objective function. 
Considering Table 1, eight parameters should be evaluated prior to employing the [WO in real 
engineering problems. Some algorithms have more parameters to tune than others such as the 
IWO. Different values were examined to determine initial population, maximum number of 
plants population, minimum number of seeds, maximum number of seeds, nonlinear modulation 
index, initial value of standard deviation, final value of standard deviation, and maximum 


number of iterations, but the most appropriate ones in solving the concrete gravity dam problem 
are shown in Table 1. 


In Fig. 2, the convergence of the objective function for M1 using the [WO is shown. As it is 
obvious, the objective function of shape optimization problem converged in 50 iterations. 


585 


ooo Average =— = — Minimum 


Maximum 


Objective Function (m2) 


Iteration 
Fig. 2. Convergence for IWO. 


Optimization results of two studied models and algorithms in previous works 1.e., Differential 
Evolution (DE), Charged System Search (CSS), Colliding Bodies Optimization (CBO), and 
Enhanced Colliding Bodies Optimization (ECBO) are shown in Table 2. In addition, Tilari Dam 
parameters are shown in aforementioned table. The cross-sectional area of this dam, which was 
constructed in India, is 709.493 (m2). As dams are normally constructed in wide valleys, small 
changes in their cross-sectional area lead to high-cost saving. M1 has the same upper and lower 
bounds as studies of [13] and [16]. According to the results, TWO with same conditions of DE, 
CSS, CBO, and ECBO could find the same objective function of them. In other words, IWO 
(M1), DE, CSS, CBO, and ECBO were succeeded in reducing total cross sectional area on Tilari 
Dam more than 20% 1.e., decrease from 709.493 (m2) to 564.496. The upstream and downstream 
slope faces and parameter x/ in these optimal models were 0.1, 0.6, and 28.96, respectively. 
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Table 1 

Parameters of WO. 
Parameter Symbol Value 
Number of initial population No 10 
Maximum number of plant population Pinay 100 
Minimum number of seeds Sinin 2 
Maximum number of seeds Smax 5 
Nonlinear modulation index n 3 
Initial value of standard deviation O vil 1 
Final value of standard deviation Final 0.001 
Maximum number of iterations ibs a 50 


In M2, gravity dam model had vertical upstream face. This situation can reduce dam's mass and 
water weight in upstream (resisting moments), and cross-sectional area. The results showed the 
gravity dam with perpendicular upstream face could lead to a far more economical design. Total 
cross-sectional area in M2 was calculated 522.56 (m2), which had about 26% reduction in 
comparison with prototype model (Tilari Dam). In M2 the calculated values of the parameters a, 
and ah were more than other optimal models. It is worth mentioning that these two parameters 
are chosen based on seismicity of dam’s zone. This issue was mentioned in [13], too. In fact, the 
more seismicity in dam’s site causes the more increase in the value of ah. In some studies is 
proposed to choose parameter a, value of 1/2 or 2/3 of ay. Generally, the magnitude of an 
earthquake depends on various parameters such as dam’s weight and type, dam's material 
behavior, and earthquake magnitude. Stability, stress, and geometry constraints (Eq. (31)) were 
applied in the current problem to ensure real-dam-design conditions. These constraints were 
between desired limits, which are shown in Eqs. (30) and (31). 


Table 2 
Parameters of prototype and optimal models. 
Algorithm 
Design variable Tilari Dam ig ng or rr 
IWO (M1) IWO(M2) DE CSS CBO ECBO 
i 0.1 0.1 0 0.1 0.1 0.1 0.1 
Ae 0.85 0.6 0.6 0.6 0.6 0.6 0.6 
x1 (m) 30.95 28.96 ; 28.96 28.96 28.96 28.96 
ay - 0.05 0.2 0.053 0.0589 0.0502 0.05 
a, - 0.05 0.1491 0.064 0.0558 0.0514 0.05 


Cross-sectional 


area (m) 709.493 564.49583 522.56175 564.496 564.49583 564.49583 564.49583 
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4. Conclusions 


Without any doubt, optimization techniques could reduce the construction time and costs. 
Among all infrastructures, gravity dams consume the sheer volume of materials. Accordingly, 
small changes in their geometries cause a major variation in the construction time and costs. In 
the current study, a nature-inspired algorithm, namely invasive weed optimization (IWO) was 
employed to optimize the shape of a concrete gravity dam. A real benchmark design problem 
(i.e., Tilari Dam, which is built in the India, with 709.493 m2 cross-sectional area) is used as a 
case study. The current framework can also be used in the future designs of other gravity dams. 
The performance of [WO was also compared to four algorithms in the literature: differential 
evolution (DE), charged system search (CSS), colliding bodies optimization (CBO), and 
enhanced colliding bodies optimization (ECBO). Various vertical and horizontal loads (i.e., 
water, seismic, wave, uplift, silt and so on) affect design of dams so the programming model 
should contain all of them. Two models were presented and their results were compared to four 
aforementioned algorithms. First model (M1) had the same conditions as the previous works. 
While the [WO had the same result compared to DE, CSS, CBO, and ECBO, convergence graph 
showed that the [WO converges faster than them. M1 and those four models in the literature 
reduced cross-sectional area of Tilari Dam approximately 20 percent. It is worth mentioning that 
in the gravity dams design, like any real design problem, there are different methods and codes 
and their consequent certain coefficients and considerations, which must be considered as well in 
order to compare the results obtained by different works. In addition, model 2 (M2) was 
proposed to evaluate a different condition i.e., a concrete gravity dam with perpendicular 
upstream face. According to the results, the cross-sectional area which was optimized with this 
assumption needed 26 percent less amount of concrete than real design. In summary, IWO could 
design a cross section for the gravity dam with a better vertical upstream face which endures all 
vertical and horizontal loads and had less concrete volume. 


Funding 


This research received no external funding. 


Conflicts of interest 


The authors declare no conflict of interest. 


References 


[1] Moffett M, Fazio MW, Wodehouse L. A world history of architecture. Laurence King Publishing; 
2003. 


[2] Chen S-H. Hydraulic Structures. Berlin, Heidelberg: Springer Berlin Heidelberg; 2015. 
doi:10.1007/978-3-662-47331-3. 


[3] | Ferdowsi A, Farzin S, Mousavi S-F, Karami H. Hybrid bat & particle swarm algorithm for 
optimization of labyrinth spillway based on half & quarter round crest shapes. Flow Meas Instrum 
2019;66:209-17. doi: 10.1016/j.flowmeasinst.2019.03.003. 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


[10] 


[11] 


[12] 


[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


A. Ferdowsi et al./ Journal of Soft Computing in Civil Engineering 4-3 (2020) 65-78 77 


Hojjati SH, Salehi Neyshabouri SAA. The objective design of triangular bucket for dam’s spillway 
using Non-dominated Sorting Genetic Algorithm II: NSGA-II. Sci Iran 2017;24:19-27. 
doi:10.24200/sci.2017.2373. 


A F, S-F M, S F, H K, M E, M A. Optimal Design and Hydraulic Performance Improvement of 
Labyrinth Spillway using Cuckoo Search Algorithm. Iran J Irrig Drain 2019;5:1086—97. 


Hassanvand MR, Karami H, Mousavi S-F. Use of multi-criteria decision-making for selecting 
spillway type and optimizing dimensions by applying the harmony search algorithm: Qeshlagh 
Dam Case Study. Lakes Reserv Res Manag 2019;24:66—75. doi:10.1111/lre.12250. 


Norouzi R, Daneshfaraz R, Ghaderi A. Investigation of discharge coefficient of trapezoidal 
labyrinth weirs using artificial neural networks and support vector machines. Appl Water Sci 
2019;9:148. doi: 10.1007/s13201-019-1026-5. 


Ehteram M, Karami H, Mousavi S-F, Farzin S, Kisi O. Evaluation of contemporary evolutionary 
algorithms for optimization in reservoir operation and water supply. J Water Supply Res Technol - 
Aqua 2018;67:54—67. doi: 10.2166/aqua.2017.109. 


Karami H, Mousavi SF, Farzin S, Ehteram M, Singh VP, Kisi O. Improved Krill Algorithm for 
Reservoir Operation. Water Resour Manag 2018;32:3353—72. doi:10.1007/s11269-018-1995-4. 


Choopan Y, Emami S. Optimal Operation of Dam Reservoir Using Gray Wolf Optimizer 
Algorithm (Case Study: Urmia Shaharchay Dam in Iran). J Soft Comput Civ Eng 2019;3:47-61. 


Rezaeeian A, Davoodi M, Jafari MK. Determination of Optimum Cross Section of Earth Dams 
Using Ant Colony Optimization Algorithm. Sci Iran 2018:0—0. doi:10.24200/sci.2018.21078. 


Ghazvinian H, Karami H, Farzin S, Mousavi SF. Effect of MDF-Cover for Water Reservoir 
Evaporation Reduction, Experimental, and Soft Computing Approaches. J Soft Comput Civ Eng 
2020;4:98—110. doi:10.22115/scce.2020.213617.1156. 


A K, P Z. Stability based optimum design of concrete gravity dam using CSS, CBO and ECBO 
algorithms. Int J Optim Civ Eng 2015;5:419-31. 


F. S. Design of gravity dam by genetic algorithms. Int J Civ Environ Eng 2011;3:187—92. 


Khatibinia M, Khosravi S. A hybrid approach based on an improved gravitational search algorithm 
and orthogonal crossover for optimal shape design of concrete gravity dams. Appl Soft Comput 
2014; 16:223—33. doi:10.1016/j.asoc.2013.12.008. 


R D, C. S. Optimal design of gravity dam using differential evolution algorithm. Int J Optim Civ 
Eng 2015;5:255—66. 

CHITI H, KHATIBINIA M, AKBARPOUR A, NASERI HR. Reliability—based design 
optimization of concrete gravity dams using subset simulation. Int J Optim Civ Eng 2016;6:329— 
48. 

Khatibinia M, Chiti H, Akbarpour A, Naseri HR. Shape optimization of concrete gravity dams 
considering dam—water—foundation interaction and nonlinear effects. [UST 2016;6:115—34. 
Memarian T, Shahbazi Y. Integrated Metaheuristic Differential Evolution Optimization Algorithm 
and Pseudo Static Analysis of Concrete Gravity Dam. Civ Eng J 2017;3:617—25. 

Zhang M, Li M, Shen Y, Zhang J. Isogeometric shape optimization of high RCC gravity dams with 
functionally graded partition structure considering hydraulic fracturing. Eng Struct 2019;179:341— 
52. doi:10.1016/j.engstruct.2018.11.005. 


Chassiakos AP, Rempis G. Evolutionary algorithm performance evaluation in project time-cost 
optimization. J Soft Comput Civ Eng 2019;3:16—29. 


78 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


[35] 


[36] 


[37] 


A. Ferdowsi et al./ Journal of Soft Computing in Civil Engineering 4-3 (2020) 65-78 


Javidrad F, Nazari M, Javidrad HR. An Innovative Optimized Design for Laminated Composites in 
terms of a Proposed Bi-Objective Technique. J Soft Comput Civ Eng 2020;4:1—28. 


Aiyelokun O, Ojelabi A, Agbede O. Performance Evaluation of Machine Learning Models in 
Predicting Dry and Wet Climatic Phases. Soft Comput Civ Eng 2020;4:29-48. 


Kalman Sipos T, Parsa P. Empirical Formulation of Ferrocement Members Moment Capacity 
Using Artificial Neural Networks. J Soft Comput Civ Eng  2020;4:111-26. 
doi:10.22115/scce.2020.221268.1181. 


Wolpert DH, Macready WG. No free lunch theorems for optimization. IEEE Trans Evol Comput 
1997; 1:67—82. doi:10.1109/4235.585893. 


Mehrabian AR, Lucas C. A novel numerical optimization algorithm inspired from weed 
colonization. Ecol Inform 2006;1:355—-66. doi:10.1016/j.ecoinf.2006.07.003. 


Asgari H-R, Bozorg Haddad O, Pazoki M, Lodiciga HA. Weed Optimization Algorithm for 
Optimal Reservoir Operation. J Irrig Drain Eng 2016;142:04015055. doi:10.1061/(ASCE)IR.1943- 
4774.0000963. 


Azizipour M, Ghalenoei V, Afshar MH, Solis SS. Optimal Operation of Hydropower Reservoir 
Systems Using Weed Optimization Algorithm. Water Resour Manag 2016;30:3995—4009. 
doi:10.1007/s11269-016-1407-6. 


Ehteram M, P. Singh V, Karami H, Hosseini K, Dianatikhah M, Hossain M, et al. Irrigation 
Management Based on Reservoir Operation with an Improved Weed Algorithm. Water 
2018;10:1267. doi:10.3390/w 10091267. 


A. F. Optimization of labyrinth spillway geometry using metaheuristic algorithms under climate 
change. M.Sc. thesis. Semnan University, Semnan, Iran. 2019. 


R S, K. P. Differential evolution - A simple and efficient adaptive scheme for global optimization 
over continuous spaces. Technical report TR-95-012. Int Comput Sci Berkeley, Calif 1995. 


Kaveh A, Talatahari S. A novel heuristic optimization method: charged system search. Acta Mech 
2010;213:267—89. doi:10.1007/s00707-009-0270-4. 


Kaveh A, Mahdavi VR. Optimal Design of Truss Structures with Continuous Variables Using 
Colliding Bodies Optimization. Colliding Bodies Optim, Cham: Springer International Publishing; 
2015, p. 39-86. doi: 10.1007/978-3-319-19659-6_3. 


Kaveh A, Ilchi Ghazaan M. Enhanced colliding bodies optimization for design problems with 
continuous and discrete variables. Adv Eng Softw 2014;77:66—75. 
doi:10.1016/j.advengsoft.2014.08.003. 


Yang X-S, editor. Nature-Inspired Algorithms and Applied Optimization. vol. 744. Cham: Springer 
International Publishing; 2018. doi:10.1007/978-3-3 19-67669-2. 


Ehteram M, Singh VP, Ferdowsi A, Mousavi SF, Farzin S, Karami H, et al. An improved model 
based on the support vector machine and cuckoo algorithm for simulating reference 
evapotranspiration. PLoS One 2019;14:e0217499. doi:10.1371/journal.pone.0217499. 


Ferdowsi A, Mousavi S-F, Farzin S, Karami H. Optimization of dam’s spillway design under 
climate change conditions. J Hydroinform 2020;22:916—36. doi:10.2166/hydro.2020.019. 


